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ABSTRACT 

We present an analysis of high resolution HI data of the dwarf galaxy DDO 47, aimed at testing the 
hypothesis that dark halo triaxiality might induce non-circular motions resulting in rotation curves 
best fitted by cored halos, even if the dark matter halo is intrinsically cuspy. This hypothesis could be 
invoked in order to reconcile the predictions of the standard ACDM theory with the rotation curves 
of disk galaxies. DDO 47 is an ideal case to test this hypothesis because it has a very regular velocity 
field, its rotation curve is best fitted by a cored halo and a NFW halo is inconsistent with the data. We 
analysed the velocity field through the higher-order harmonic terms in order to search for kinematical 
signatures of alleged non-circular motions needed to "hide" a cusp: the result is that globally non- 
circular motions are at a level of 2-3 km s" 1 , and they are more likely to be associated to the presence 
of some spiral structure than to a global elongated potential (e.g. a triaxial halo). These non-circular 
motions are far from being sufficient to account for the discrepancy with the ACDM predictions. We 
therefore conclude that the dark matter halo around the dwarf galaxy DDO 47 is truly cored and that 
a cusp cannot be hidden by non-circular motions. 

Subject headings: galaxies: dwarf - galaxies: individual (DDO 47) - galaxies: kinematics and dynamics 
- dark matter. 



1. INTRODUCTION 

A fundamental prediction of the cosmological (A) Cold 
Dark Matter (ACDM) theory is that virialized dark mat- 
ter halos should have universal density profiles. An em- 
pirical fit formula for the spherically averaged density 
distribution p(r) had been proposed by Navarro, Frenk 
& White (1996, 1997, NFW) 



p(r) = 



{r/r s )(l + r/r s ) 2 



(1) 



However, as we can typically trace rotation curves ob- 
servationally in the regime of 0.1r s < r < 2r s , the NFW 
profiles still provide a good description to the dark mat- 
ter halos in galaxies. 

The CDM predictions have been confronted with ob- 
servations, investigating the rotation curves of disk 
galaxies. Numerous mass models of spiral, dwarf and 
LSB galaxies have shown that in a large number of cases 
the inferred cores of dark matter halos are much shal- 
lower than expected from the NFW-fit (Flores & Pri- 
mack 199 4, Moore 1994, Burkert 1995, Salucci & Burkert 



where p s and r s are two in principle free but for a given 
cosmological model strongly correlated parameters. This 
density distribution converges to power-law functions at 
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et al. (1999) later on argued that the inner cusp may even 
be steeper with p ~ r~ 15 . 

The power-law behaviour has been a matter of strong 
debate. Recent simulations, e.g. by Power et al. (2003) 
and Reed et al. (2005) demonstrate that the CDM halo 
density profiles are better represented by a function with 
a continuously varying slope and no convergence to a 
power-law slope at small radii. A revised fitting formula 
was derived by Navarro et al. (2004), where the slope 
at the innermost radius is about -1.2 for spiral galaxies 
and -1.35 for dwarfs. Merritt et al. (2005) showed that 
the Sersic (1968) law fits high-resolution ACDM halos 
with a wide range of masses with a Sersic index of n w 
2.4 — 3 which is usually used as a good fit of the surface 
brightness profiles of the centres of early-type galaxies. 
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(1 +r/r c )(l + r/r c ) 2 



(2) 



fits the data often very well, indicating that, in con- 
trast to cosmological predictions, the maximum density 
of dark halos is finite. Donato et al. (2004) even found 
a strong correlation between the core radius r c and the 
dynamical mass indicating again that the observed cored 
dark matter halos represent a 1-parameter family. 

This contradiction between theoretical predictions and 
observations has stimulated a lot of discussions as it has 
the potential to provide interesting new insights into the 
nature of dark matter and its possible interaction with 
visible matter (for reviews see Primack 2004 or Ostriker 
& Steinhardt 2003). In this letter we will explore the 
possibility that solid-body rotation curves which are usu- 
ally interpreted as a signature of constant density cores 
would arise naturally also in NFW halos as a result of 
a triaxial dark matter mass distribution that hides the 
steep central density cusp inducing deviations from ax- 
ial symmetry in the disk (e.g. Hayashi et al. 2004); see 
however Gnedin et al. 2004 and Kazantzidis et al. 2004 
for the effect of baryons reducing halo triaxiality. The 
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7 h 39 m 20 5 39 m 10 s 39 m 00 s 38 m 50 s 

Right Ascension (B1950) 

Fig. 1.— Velocity field of DDO 47 (contours) from the 15 " 
resolution cube and its total HI map (grey scale). Contours are 
spaced by 10 km s — 1 . 

deviations from axial symmetry in the gas kinematics 
can be tested observationally with galaxies showing two- 
dimensional velocity fields that are clearly inconsistent 
with a NFW profile. 

2. THE TARGET: DDO 47 

An ideal case is the nearby dwarf galaxy DDO 47 for 
which HI observations (Walter & Brinks 2001) with ad- 
equate resolution and sensitivity exist; these data show 
that the kinematics of the gas is very regular, with a 
well-behaved velocity field (Fig. 1). 

DDO 47 is one of the most clear cases with kinematical 
properties that are inconsistent with NFW. A first rota- 
tion curve of DDO 47 was presented in Walter & Brinks 
(200f ) and a more detailed decomposition, performed by 
Salucci, Walter & Borriello (2003), shows that a NFW 
halo is inconsistent with the data. They also show that 
the halo of DDO 47 has a central constant-density core: 
the rotation curve was successfully fitted with a Burkcrt 
halo with r c =7 kpc and po = 1.4 x 10~ 24 g cm~ 3 (as- 
suming a distance of 4 Mpc, and 1 arcmin corresponds to 
1.16 kpc); any other cored halo (e.g. pseudo- isothermal) 
gives equally good fits to the rotation curve. We show 
the Burkert fit in Fig. 2, where the rotation curve was 
derived by leaving the position angle as a free parame- 
ter, in order to better trace the warp. The error bars are 
derived from the difference between the two sides, also 
considering a minimum error of half the velocity resolu- 
tion, corrected by the inclination. Also shown in Fig. 2 
is the NFW fit; differently from Salucci, Walter & Bor- 
riello (2003) we do not fix the virial mass M vir by match- 
ing the last velocity point, but we leave it instead as a 
free parameter, the concentration parameter c scaling as 
c = 21( 10 ^J -°- 13 (Wechsler et al. 2002). 

Here we subject the same data to a more thorough 
analysis, carried out in the light of the above discussion. 
The spatial resolution of the datacubc (~ 300 pc) is sim- 
ilar to that of the dwarf galaxy simulations performed by 
Navarro et al. (2004). 




100 150 200 

r (arcsec) 

Fig. 2. — Mass models for DDO 47: the solid line is the best- 
fitting model with the Burkert halo, and the dotted line represents 
the best-fitting NFW model, with a virial mass M v i r = 2.4 X 
1O 1O M (see text). 
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Fig. 3. — Tiltcd-ring fit of the velocity field with the inclination 
as a free parameter. From top to bottom: inclination, position 
angle (fixed in this case, see text) and velocity as a function of 
radius. The uncertainties shown are formal errors. 

We have made a tilted-ring analysis of the velocity field 
with the inclination and the centre fixed at the same val- 
ues as used by Walter & Brinks (2001). As noticed by 
these authors, leaving the inclination i as a free param- 
eter does not give a very stable result. In order for a 
possible radial change of i to account for the discrep- 
ancy with the ACDM predictions, i should increase as a 
function of radius. Instead, as is shown in Fig. 3, the 
tilted-ring fit with the inclination free (with 35° as initial 
guess) suggests that the radial change of the inclination, 
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Fig. 4. — Harmonic decomposition of the velocity field, with 
terms up to the order 3. Notice that here the coefficients of cq. 
3 with j > 1 have been divided by sin i in order to have physical 
velocities instead of line of sight velocities. The uncertainties shown 
are formal errors. 

if any, goes in the opposite direction (i.e. i seems to 
get larger and not smaller for small radii). Note that in 
this figure we have fixed the position angle at an average 
value (separately for each radius) between fits made with 
different values of the inclination, between 5° and 65°. 
Therefore, considering that there is no evidence for an 
increase of the inclination, we conclude that fixing the 
inclination at 35° is a reasonable assumption for the aim 
of this paper. We have also found that a slight change in 
the systemic velocity (273 km s^ 1 instead of 272 km s _1 ) 
decreases the differences between the two sides. 

3. HARMONIC DECOMPOSITION, NON-CIRCULAR 
MOTIONS AND HALO TRIAXIALITY 

In order to better investigate the kinematic properties 
of DDO 47, we have fitted the velocity field in terms of 
harmonic coefficients, where the observed velocity along 
the line of sight Vi os is given by: 

n 

Vios = c +^2[cjCOs(ji>) + SjSin(jV)] (3) 
i=i 

where ip is the azimuthal angle defined as the angle 
from the major axis of the receding side in a counterclock- 
wise direction. In this analysis we only considered terms 
up to j = 3. In the case of a rotating disk plus axisym- 
metric radial motion (i.e. the traditional tilted-ring fit 
of the velocity field), only terms up to j = 1 are present 
and Co = V sys (the systemic velocity), ci = Kot sin i 
(the rotation velocity, and i is the inclination angle), 
and si = Kad sm i (the radial velocity). More detailed 
descriptions of this method and its applications can be 
found in Schoenmakers, Franx and de Zeeuw (1997), 
Schoenmakers (1999) and Wong, Blitz & Bosma (2004). 

The results of the harmonic decomposition of the ve- 
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Fig. 5. — Residuals of the tilted-ring fit. Contours are -15, -10, 
-5, 5, 10, and 15 km s _1 . 
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Fig. 6. — Observed rotation curves of the dwarf galaxy DDO 47 
with the position angle fixed, along the major axis (filled squares) 
and along the minor axis (open squares). The stellar exponential 
scale length is 0.5 kpc. 

locity field are shown in Fig. 4, where the parameters 
were fitted with a fixed inclination (35°) and the posi- 
tion angle left as a free parameter. First of all it can 
be noticed that the s\ term is negligible and that the S3 
term is larger but still has a small amplitude, at most 3 
km s _1 . The fact that S3/S1 >> 1 is to be expected, as 
DDO 47 has a low inclination and this fraction is a strong 
function of inclination (Schoenmakers 1999). As shown 
by Schoenmakers et al. (1997), a global elongation of the 
potential would cause the S3 term to be approximately 
constant and non-zero, while wiggles could be explained 
by the presence of spiral structure. In DDO 47 the latter 
case applies, we find no evidence for a global elongation 
of the potential and the amplitude of the oscillations is 
much smaller (a factor 5-10) than the discrepancy with 
the ACDM predictions. Fig. 4 also shows that some 
lopsidedness is present, since the j » = 2 terms are non- 
zero. In the residual velocity field (Fig. 5) we find that 
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the regions with residuals larger than 5 km _1 (9 km _1 
deprojected) are not a common feature. The innermost 
ones are likely to be associated to either star-forming re- 
gions or to the possible bar, while the outermost ones are 
probably due to asymmetry in the outer warp. 

As can be seen in Walter & Brinks (2001), the ma- 
jor axis of the stellar body seems to be a bit misaligned 
with the HI disk, and low-level non-circular motions as- 
sociated to the optical part of the galaxy can be seen in 
the velocity field. Unfortunately due to the large number 
of stars in the field a more thorough analysis of the pos- 
sible bar structure in the optical image is not possible. 
If the S3 ~ 3 km s _1 region around 60-70" were asso- 
ciated to a barred potential, its amplitude would any- 
way be much lower than needed to reconcile the rotation 
curve with the ACDM predictions. We have discussed 
so far two of the three effects that according to Rhee 
et al. (2004) could bias the rotation curve (the inclina- 
tion and the bar). Concerning the third effect, i.e. the 
influence of a possible bulge, these authors claim that 
the random motions of the bulge material could bias the 
observed rotation velocities towards the systemic veloc- 
ity. However in DDO 47 visual inspection of the velocity 
profiles shows that they are symmetric. This was further 
assessed by building a velocity field where the velocity 
at any point is given by the peak of the profile: if the 
profiles were strongly asymmetric, the difference between 
this peak velocity field and a first-moment velocity field 
would be large and the separation between approaching 
and receding side should be evident. Instead, we see 
small deviations and no obvious separation between the 
two sides. 

We can therefore conclude that the non-circular mo- 
tions in DDO 47 are at a level of at most 3 km s _1 and 
that they are not associated with a global elongation of 
the potential. As a particular case in Fig. 6 we show the 
minor axis kinematics of DDO 47: the velocity shows 
very small deviations from zero; the analysis performed 



earlier gives a more general result, making use of the 
whole velocity field to derive the amount of non-circular 
motions. We conclude that the global discrepancy be- 
tween observations and ACDM predictions shown in Fig. 
2 cannot be explained by non-circular motions induced 
by halo triaxiality, regardless of the viewing angle. 

4. CONCLUSIONS 

We have presented a dynamical analysis of HI data of 
the dwarf galaxy DDO 47 to test the hypothesis that 
the possible triaxial shape of the dark matter halo might 
induce deviations from circular motions resulting in a ro- 
tation curve best fitted by a cored halo, with the dark 
matter halo being still CDM-like (i.e., cuspy). The high- 
quality HI observations of DDO 47 make it an ideal case 
to test this hypothesis because it has a very regular veloc- 
ity field, its rotation curve is best fitted by a cored halo 
and a NFW halo is inconsistent with the data (Salucci 
et al. 2003). Even though our analysis indicates a cored 
dark matter halo profile for DDO 47, work to test the 
effect of different halo shapes and disk/halo alignments 
on the shapes of rotation curves is clearly very valuable. 

We have performed a harmonic decomposition of the 
velocity field in order to estimate the amount of non- 
circular motions. We find that they are at a level of at 
most 3 km s _1 and that they are likely to be associated 
to spiral structure rather than to a global elongation of 
the potential. In any case the observed deviations from 
circular motion are much smaller than the discrepancy 
between ACDM predictions and the observed rotation 
curve. We conclude therefore that the dark matter halo 
around the dwarf galaxy DDO 47 is truly cored and that 
a cusp cannot be hidden by non-circular motions. 
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